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1 Introduction 



Successfully describing heavy quarkonium production is a long-standing problem in QCD. 
An effective field theory, non-relativistic QCD (NRQCD) [lj[2], provides the foundation 
for much of the current theoretical work. According to NRQCD, the production of heavy 
quarkonium factorises into two steps: a heavy quark-antiquark pair is first created at 
short distances and subsequently evolves non-perturbatively into quarkonium at long 
distances. The NRQCD calculations depend on the colour-singlet (CS) and colour-octet 
(CO) matrix elements, which account for the probability of a heavy quark-antiquark pair 
in a particular colour state to evolve into a heavy quarkonium state. The CS model 
(CSM) (3||4), which provides a leading-order description of quarkonium production, was 
initially used to describe experimental data. However, it underestimates the observed cross- 
section for single J/tp production at high transverse momentum (pp) at the Tevatron [5). 
To resolve this discrepancy, the CO mechanism was introduced [6J. The corresponding 
matrix elements were determined from the high-px data, as the CO cross-section decreases 
more slowly with px than that predicted by CS. More recent higher-order calculations (7 10 



close the gap between the CS predictions and the experimental data |11 , reducing the 
need for large CO contributions. 

Studies of the production of the J/ip and T(1S), T(2S) and T(3S) mesons (indicated 
generically as T in the following) have been performed using pp collision data taken at 

y/s = 7 TeV and at = 2.76 TeV by the LHCb (T^H' ALICE @0' ATLAS HUE 
and CMS 20 22 experiments in different kinematic regions. As well as providing direct 
tests of the underlying production mechanism, these studies are crucial to estimate the 
contribution of double parton scattering to multiple quarkonium production 23 24]. 

In this paper first measurements of quarkonium production at \/s = 8 TeV are 
reported under the assumption of zero polarisation, an assumption that is discussed 
in the paper. The differential production cross-sections of prompt J/ip and T mesons, 
produced at the pp collision point either directly or via feed-down from higher mass 
charmonium or bottomonium states, are presented in the range of rapidity 2.0 < y < 4.5 
and pt < 14 GeV/c (J/ip) or px < 15 GeV/c (T). The fraction of J/ip mesons from 
6-hadron decays, abbreviated as u J/ip from 6" in the following, is also measured in the 
same fiducial region. 



2 The LHCb detector and data set 



The LHCb detector 25 is a single-arm forward spectrometer covering the pseudorapidity 
range 2 < rj < 5, designed for the study of particles containing b or c quarks. The detector 
includes a high precision tracking system consisting of a silicon-strip vertex detector 
surrounding the pp interaction region, a large-area silicon-strip detector located upstream 
of a dipole magnet with a bending power of about 4 Tm, and three stations of silicon-strip 
detectors and straw drift tubes placed downstream. The combined tracking system has 
a momentum resolution Ap/p that varies from 0.4% at 5GeV/c to 0.6% at lOOGeV/c, 
and an impact parameter resolution of 20 \mi for tracks with high px- Charged hadrons 
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are identified using two ring-imaging Cherenkov detectors |26|. Photon, electron and 
hadron candidates are identified by a calorimeter system consisting of scintillating-pad and 
preshower detectors, an electromagnetic calorimeter and a hadronic calorimeter. Muons 
are identified by a system composed of alternating layers of iron and multiwire proportional 
chambers, with the exception of the centre of the first station, which uses triple-GEM 
detectors. 

The data sample used in this analysis was collected during the first part of the data 
taking period at \fs = 8 TeV in April 2012. During this period the average number 
of interactions per crossing varied. The T meson analysis is based on a data sample, 
corresponding to an integrated luminosity of about 51 pb _1 of pp interactions, collected 
with an average of 1.3 visible interactions per crossing. The analysis for the more abundant 
J/ip mesons is based on data, corresponding to an integrated luminosity of about 18 pb _1 , 
collected with an average of 1.0 visible interactions per crossing. The trigger 27 consists of 
a hardware stage, based on information from the calorimeter and muon systems, followed 
by a software stage, which applies a full event reconstruction. At the hardware stage, 
events are selected requiring dimuon candidates with a product of their pp larger than 
1.68 ( GeV/c) 2 . In the subsequent software trigger, two well reconstructed tracks are 
required to have hits in the muon system, a higher than 500 MeV/c, p higher than 
6 GeV/c and to form a common vertex. Only events with a dimuon candidate with an 

o2 



invariant mass m w within 120 MeV/c of the known J/ip meson mass 28] or larger than 



4.7 GeV/c are retained for further analysis. 



3 Selection and cross-section determination 



The selection is based on the criteria described in Refs. 12, 13 and is summarised in 
Table [T] It starts by combining oppositely-charged particles, identified as muons, with 
a track pp larger than 700 (1000) MeV/c 2 for the J/i/j(T) meson. Good track quality is 
ensured by requiring a \ 2 per degree of freedom, x 2 /ndf, less than 4 for the track fit. 
Duplicate particles created by the reconstruction are suppressed to the level of 0.5 x 10~ 3 
using the Kullback-Leibler (KL) distance variable 29 . To ensure good quality vertex 



'0. 



reconstruction, the x probability of the dimuon vertex is required to be larger than 0.5 ' 
In addition, the primary vertex (PV) associated to the dimuon candidate is required to be 
within the luminous region, defined as |xpv| < 1 nim , |?/pv| < 1 nim and |zpv| < 150 mm. 

In the J/0 analysis additional criteria are applied to the vertex quality. The uncertainty 
on the pseudo decay time t z , defined in Eq. |2j is required to be less than 0.3 ps, as estimated 
by the propagation of the uncertainties given by the track reconstruction. 

The simulation samples are based on the Pythia 6.4 generator 



the parameters detailed in Ref. 31 . The EvtGen package 32 is used to generate hadron 



30 configured with 



decays. The interaction of the generated particles with the detector and its response are 

Radiative corrections 



implemented using the Geant4 toolkit 33 as described in Ref. 34 



to the decay of the vector meson to dimuons are generated with the Photos package (35] . 
The differential cross-section for the production of a vector meson V in a bin of (pp, y), 



2 



Table 1: Selection criteria for the J/ip and T meson analyses. Criteria common to both analyses 
are displayed between the two columns. 



Variable Value (J/0) Value (T) 

Track p T ( MeV/c) > 700 > 1000 

Track xV n df < 4 

KL distance > 5000 

Vertex \ 2 probability > 0.5% 

t z uncertainty (ps) < 0.3 

Mass window m w ( MeV/c 2 ) |m MM - M(J/ip ) | < 120 8500 < m m < 11500 



where V stands for a J/ip or T meson, decaying into a muon pair, is 

d 2 cr = NjV^n+fj,-) 

dydp T [PP ^ j £xe tot xB(F^//+/i-)xAyxAp T ' U 

where N (V — > is the number of observed V — > candidates, e to t the total 

detection efficiency in the given bin, C is the integrated luminosity, B (V — > is the 

branching fraction of the V — > fi + ^~ decay and Ay = 0.5 and Apx = 1 GeV/c are the 
rapidity and pr bin sizes, respectively. In the case of the J/ip — > decay the branching 



fraction is well known, B(J/ijj — > ) = (5.94 ± 0.06) x 10 28 , and therefore it 



is chosen to quote an absolute cross-section. On the other hand, the dimuon branching 



fractions of the T mesons are known less precisely [28], and therefore, as in Ref. 13 , the 
product of the cross-section times the dimuon branching fraction is given. 

The total efficiency e to t is the product of the geometric acceptance, the reconstruction 
and selection efficiency and the trigger efficiency. All efficiency terms are evaluated using 
simulated samples and validated with data-driven techniques in each (pt,v) bin. 

The procedure to measure the integrated luminosity is described in Ref. |36|. For this 
analysis a van der Meer scan [37] was performed in April 2012, resulting in a measurement 
of the integrated luminosity of 18.4 ± 0.9 pb~ for the J/ip and 50.6 ± 2.5 pb~ for the T 
samples. 



4 J/jJj meson signal 

As in the previous studies, prompt J/ip mesons are distinguished from J/ip from b by 
means of the pseudo decay time variable defined as 

t = (zjm ~ ^pv) x Mj/^ 

Pz 
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Figure 1: Projections of the fit result for a selected bin in px and y for (a) the J/0 invariant 
dimuon mass and (b) t z . For the former, the total fitted function is shown (blue solid line) 
together with the signal distribution (red dotted line). In the t z projection the total fitted 
function is shown together with the J/0 from b component, the prompt signal, the background 
and the tail component due to the association of a J/0 candidate with a wrong PV. 



where z j/^ and zpy are the positions along the beam axis z of the J/0 decay vertex and of 
the primary vertex refitted after removing the decay muons of the J/0 candidate; p z is the 
measured J/0 momentum in the beam direction and Mj/^ is the known J/0 mass 28 . 



The yields of both prompt J/0 mesons and J/0 from b are determined from a two- 
dimensional fit in each (pt,v) bin to the distributions of invariant mass and pseudo decay 



time of the signal candidates, following the approach described in Ref. 12 . The mass 
distribution is modelled with a Crystal Ball function [38] for the signal and an exponential 
function for the combinatorial background. 

The signal pseudo decay time distribution is described by a ^-function at t z = for the 
prompt J/0 component together with an exponential decay function for the J/0 from b 
component. The shape of the tail arising from the association of a J/0 meson candidate 
with a wrong primary vertex is derived from the data by combining a J/0 meson from a 
given event with the primary vertex of the following event in the sample. The prompt 
component of the signal function and that from b hadron decays are convolved with a 
resolution function modelled by the sum of two Gaussian functions. The background 
distribution is parameterised with an empirical function based on the shape of the t z 
distribution observed in the J/0 mass sidebands. It is built as the sum of a 5- function and 
five exponential components, three for positive t z and two for negative t z , the negative and 
positive exponential functions with the largest lifetime having their lifetimes tl fixed to 
the same value. This function is convolved with a resolution function modelled by the sum 
of two Gaussian functions. All parameters of the background component are determined 
independently in each (pr,y) bin from the distribution of the pseudo decay time and are 
fixed in the final fit. The total fit function is the sum of the products of the mass and t z 
fit functions for the signal and background. Figure [T] shows the fit projections in mass and 
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Figure 2: Invariant mass distribution of the selected T — > candidates in the range 

Pt < 15 GeV/c and 2.0 < y < 4.5. The three peaks correspond to the T(1S), T(2S) and T(3S) 
meson signals (from left to right). The superimposed curve and the signal yields (dotted) are the 
result of the fit described in the text. 



t z for one specific bin (3 < pt < 4 GeV/c, 2.5 < y < 3.0) with the fit result superimposed. 
Summing over all bins, a total signal yield of 2.6 million J/ip events is obtained. 

5 T meson signal 

The T meson signal yields are determined from a fit to the reconstructed dimuon invariant 
mass of the selected candidates with 8.5 < m MAt < 11.5 GeV/c 2 . The distribution is 
described by the sum of three Crystal Ball functions, one for each of the T(1S), T(2S) 
and T(3S) signals, and an exponential function for the combinatorial background. The 
parameters a and n of the Crystal Ball function describing the radiative tail are fixed to 
the values of a = 2 and n = 1 based on simulation studies. The width of the Crystal Ball 
function describing the T(1S) meson is allowed to vary, while the widths of the T(2S) 
and T(3S) mesons are constrained to the value of the width of the T(1S) signal, scaled 
by the ratio of the masses of the T(2S) and T(3S) to the T(1S) meson. The peak values 
of the T(1S), Y(2S) and T(3S) mass distributions are allowed to vary in the fit and are 
consistent with the known values [28] . Figure [2] shows the results of the fit performed over 
the entire range in and y. The obtained signal yields are 43 785 ± 254, 10 976 ± 155 and 
5325 ± 122 for the T(1S), T(2S) and T(3S) mesons, respectively, with a mass resolution 
of the T(1S) resonance of 43MeV/c 2 . The fit is repeated independently for each of the 
bins in p^ and y. When fitting the individual bins, the masses are fixed to the values 
obtained when fitting the full range, while the mass resolution for the T(1S) candidates is 
parameterised with a linear function of px and y that reproduces the behaviour observed 
in data. Bins with too few entries are excluded from the analysis. 
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Table 2: Relative systematic uncertainties (in %) on the J/tp and T cross-section results and on 
the fraction of J/0 from b. 



Correlated between bins 



Mass fits 


0.7 to 2.2 


Radiative tail 


1.0 


Muon identification 


1.3 


Tracking efficiency 


0.9 


Vertexing 


1.0 


Trigger 


4.0 


Luminosity 


5.0 


B(J/ip ->• n+fi~) 


1.0 


Uncorrelated between bins 


Production model 


1.0 to 6.0 


t z fit, for J/ip from b 


1.0 to 12.0 



6 Systematic uncertainties 

Previous studies [l~2~|[l~3~] have shown that the total efficiency depends on the initial 
polarisation state of the vector meson. The J/0 polarisation has been measured at 

and ALICE (40] collaborations, in a kinematic range similar 
to that used in this analysis, and the T polarisation has been measured by CMS [4l] 
at large and central rapidity. They were both found to be small. Therefore, in this 
paper results are quoted under the assumption of zero polarisation and no corresponding 
systematic uncertainty is assigned on the cross-section for this effect. All other systematic 
uncertainties are summarised in Table [U 

Uncertainties related to the mass model describing the shape of the dimuon mass 
distribution are estimated by fitting the invariant mass distributions for the J/0 and T 
mesons with the sum of two Crystal Ball functions. The relative difference in the number 
of signal events (0.7-2.2%) is taken as a systematic uncertainty. A fraction of events 
has a lower invariant mass because of the energy lost through bremsstrahlung. Based on 
simulation studies, about 4% of the signal events are estimated to be outside the analysis 
mass windows and are not counted as signal. The fitted signal yields are corrected for this 
effect and an uncertainty of 1% is assigned to the cross-section measurement based on a 
comparison between the radiative tail observed in data and simulation. 

The uncertainty due to the muon identification efficiency is measured on data using a 
tag-and-probe method. This method reconstructs J/0 candidates in which one muon is 
identified by the muon system ( "tag" ) and the other ( "probe" ) is identified selecting a 



= 7 TeV by the LHCb 39 
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track depositing the energy of minimum-ionising particles in the calorimeters. 

The ratio of the muon identification efficiency measured in data to that obtained in 
the simulation is convolved with the momentum distribution of muons from J/ip and T 
mesons to obtain an efficiency correction. This is found to be 0.98 ± 0.01; the uncertainty 
on the correction factor is considered as a systematic uncertainty. 

The uncertainty on the reconstruction efficiency of the muon tracks has also been 
estimated using a data-driven tag-and-probe approach based on partially reconstructed 
J/ijj decays, and it was found to be 0.9% per muon pair. 

Differences between data and simulation in the efficiency of the requirement on the 
vector meson vertex \ 2 probability lead to a further uncertainty of 1%. 

The trigger efficiency is determined using a data-driven method exploiting a sample of 



events that are still triggered when the signal candidate is removed 27 . The efficiency 



obtained with this method in each (px, y) bin is used to check the efficiencies measured in 
the simulation. The systematic uncertainty associated with the trigger efficiency is the 
difference between that measured in the data and in the simulation. As a cross-check, the 
trigger efficiency is also computed using a data sample that has not been required to pass 
any physics trigger. The results obtained with the two methods are consistent. 

The luminosity is determined with an uncertainty of 5%, dominated by differences in 
the results obtained with a van der Meer scan 37 using the core and off-core parts of the 
beam. 

The dependence of the efficiency calculation on the production model used in the 
simulation is taken into account by varying the main parameters of the Pythia 6.4 
generator related to prompt vector meson production. These parameters define the 
minimum p^ cut-offs for regularising the cross-section. This effect is evaluated in each 
(pT,y) bin and found to be at most 6%. 

Uncertainties related to the t z fitting procedure for the J/ip mesons are included by 
changing the parameterisation used to describe the signal and background. A second 
fitting method based on the sPlot technique 42 is used with the mass as the control 
variable to unfold the background and to perform an unbinned likelihood fit to the pseudo 
decay time distribution. The two approaches give consistent results and their difference is 
taken as an estimate of the systematic uncertainty. These uncertainties are evaluated in 
each (p T , y) bin and found to be a few percent. 



7 Results on J/0 meson production 

The measured double-differential production cross-sections for prompt J/0 mesons, under 
the assumption of zero polarisation, and for J/0 from b are given in bins of and V i n 
Tables [4] and [5j respectively, and are displayed in Fig. |3j The integrated cross-section for 
prompt J/0 meson production in the defined fiducial region, summing over all bins of the 
analysis, is 

a (prompt J/ip,p T < 14 GeV/c, 2.0 < y < 4.5) = 10.94 ± 0.02 ± 0.79 ]ib, 
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where the first uncertainty is statistical and the second is systematic, computed taking 
correlations into account. The integrated cross-section for the production of J/ip from b in 
the same fiducial region is 

a (J/?p from b, p T < 14 GeV/c, 2.0 < y < 4.5) = 1.28 ± 0.01 ± 0.11 ^b. 

The total bb production cross-section is computed as 

r a (J/0 from 6, p T < 14 GeV/c, 2.0 < y < 4.5) 
a( PP ttX) = a 47r mb^l/W) ' (3) 
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Figure 3: Differential production cross-section for (a) prompt J/vp mesons and (b) J/ip from b as 
a function of pr in bins of y. It is assumed that prompt J/ip mesons are produced unpolarised. 
The errors are the quadratic sums of the statistical and systematic uncertainties. 
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Table 3: Differential production cross-section da/dy in nb for prompt J/ip mesons (assumed 
unpolarised) and for J/ip from b, integrated over pp- The first uncertainty is statistical, the 
second (third) is the part of the systematic uncertainty that is uncorrelated (correlated) between 
bins. 





y 


Prompt J/ip 


J/ip from b 


2.0 


-2.5 


5140 ±26 ±49 ±368 


717±6±12±51 


2.5 


-3.0 


5066 ± 14 ± 30 ± 363 


666±3± 9 ±48 


3.0 


-3.5 


4573 ±11 ±20 ±328 


538±3± 8 ±39 


3.5 


-4.0 


3940 ±11 ±24 ±282 


388±2± 4 ±28 


4.0 


-4.5 


3153 ±12 ±16 ±226 


240±2± 3 ±17 



where the factor = 5.4 is an extrapolation factor of the cross-section from the measured 
to the full kinematic region. This factor is obtained using the simulation as described in 
Sect, g The inclusive b^J/ipX branching fraction is B(b -> J/ipX) = (1.16 ± 0.10)% (28]. 
The resulting total bb cross-section is a(pp — > bbX) = 298 ± 2 ± 36 \xb, where the first 
uncertainty is statistical and the second is systematic, which includes the uncertainty on 
B(b — > J/ipX). No systematic uncertainty has been included for the extrapolation factor 
estimated from the simulation. For comparison, the value of the extrapolation factor 



given by NLO calculations is 5.1 43 



Tableland Fig. [4] show the differential production cross-section da/dy integrated 
over prj for unpolarised prompt J/ip mesons and J/ip from b. For both components, 
the cross-section decreases significantly between the central and forward regions of the 
acceptance. 

Table [6] and Fig. [5] give the values of the fraction of J/ip from b in the different bins of 
Pt and y. The fraction of J/ip mesons from 6-hadron decays increases as a function of pr? 
and, at constant px> decreases with increasing y, as seen in the study at = 7 TeV [12|. 



8 Results on T meson production 

The double- differential production cross-sections times the dimuon branching fractions for 
the T mesons in bins of j>t and y are given in Tables [7j [8j and [9] with the assumption of no 
polarisation. The double-differential cross-sections are displayed in Fig. [6] The integrated 
cross-sections times dimuon branching fractions B lS = B(T(iS) — > fifi), with i — 1, 2, 3, in 
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the range pr < 15 GeV/c and 2.0 < y < 4.5 are measured to be 



a{pp T(1S)X) x B 1S 
a{pp ->• T(2S)X) x B 2S 
a{pp -> r(35)X) x B 3S 



3.241 ±0.018 ±0.231 nb, 
0.761 ±0.008 ±0.055 nb, 
0.369 ± 0.005 ± 0.027 nb, 



where the first uncertainty is statistical and the second systematic. The cross-section 
times dimuon branching fractions for the three T states are compared in Fig. [7] as a 
function of p^ and y. These results are used to evaluate the ratios R lS / ls of the T(2S) to 
T(1S) and T(3S) to T(1S) cross-sections times dimuon branching fractions. Most of the 
uncertainties cancel in the ratio, except those due to the size of the data sample, to the 
model dependence and to the choice of the fit function. The ratios R lS / ls as a function 
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Figure 4: Differential production cross-section as a function of y integrated over p?, for (a) 
unpolarised prompt J/ip mesons and (b) J/tp from b. The errors are the quadratic sums of the 
statistical and systematic uncertainties. 
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Figure 5: Fraction of J/ip from 6 as a function of pt, in bins of y. 
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Figure 6: Double-differential cross-sections times dimuon branching fractions as a function of 
in bins of y for (a) the T(1S), (b) T(2S) and (c) T(3S) mesons. 
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Figure 7: Differential production cross-sections for Y(1S), Y(2S) and Y(3S) mesons times 
dimuon branching fraction (a) as a function of px integrated over y, and (b) as a function of y 
integrated over pr- 



of px and V are given in Tables 10 and 11, respectively, and shown in Fig. [8j with the 
assumption of no polarisation. For this measurement the pr range has been restricted to 
pr < 14 GeV/c and the y range to 2.0 < y < 4.0 to ensure enough counts for the three T 
states in all bins. The ratios are constant as a function of y and increase as a function of 



Pr? m agreement with previous observations by LHCb 13 and as reported by ATLAS 19 
and CMS |2ll at Js = 7 TeV. 
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Figure 8: Ratio of the Y(2S) to T(1S) and T(3S) to T(1S) cross-sections times dimuon 
branching fractions (a) as a function of px integrated over y, and (b) as a function of y integrated 
over px • 
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9 Comparison with theoretical models 



The measured differential cross-sections for the production of prompt J/0 mesons as 
a function of px axe compared in Fig. [9] to three theoretical models that assume no 
polarisation. The considered models are 

- an NRQCD model at next-to-leading order (NLO). The colour-octet matrix elements 



in this case are determined from a global fit to HERA, Tevatron and LHC data 44 45 



- an NNLO* CSM (9j[l0] ; the notation NNLO* indicates that the calculation at 
next-to-next leading order is not complete and neglects part of the logarithmic terms; 

- an NLO CSM jjj with the input parameters related to the choice of scale and charm 
quark mass given in Ref. 



44 



In these comparisons it should be noted that the predictions are for direct J/0 meson 
production, whereas the experimental measurements include feed-down from higher char- 
monium states. In particular, the contribution from J/0 mesons produced in radiative 
Xc decays in the considered fiducial range was measured to be at the level of 20% at 



y/s = 7 TeV 46 . Allowing for this contribution, as was seen in the previous studies, both 
the NNLO* CSM and the NLO NRQCD models provide reasonable descriptions of the 
experimental data. In contrast, the CSM at NLO underestimates the cross-section by an 
order of magnitude. 

The results for the production of J/0 from b can be compared to calculations based 



on the FONLL formalism 43,47 . This model predicts the 6-quark production cross- 
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Figure 9: Comparison of the differential cross-section for the production of prompt J/ip meson 
(under the assumption of zero polarisation) as a function of px with direct production in an 
NLO NRQCD model [44j|45] (orange diagonal shading), an NNLO* CSM [To] (solid yellow) and 
an NLO CSM [7] (blue vertical shading). The points show the measurements reported in this 
analysis. 
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Figure 10: Differential production cross-section for J/tjj from b (a) as a function of in the 
fiducial range 2.0 < y < 4.5, and (b) as a function of y in the fiducial range pt < 14 GeV/c. The 



FONLL prediction 43 47J is shown in yellow. The points show the measurements reported in 
this analysis. 



section, and includes the fragmentation of the 6-quark into 6-hadrons and their decay 
into J/ip mesons. In Fig. 10 the data for the differential production cross-section as a 
function of and y at y/s = 8 TeV are compared to the FONLL predictions. Good 
agreement is observed. The prediction for the total cross-section in the fiducial range of 
this measurement is 1.34lo49 M-b, in good agreement with the result presented here. In 
the measurements of the cross-section for J/ip from b at 



Fig. 
and 



11 



7 



12 



s = 2.76 [14 

TeV are compared to the FONLL prediction. The behaviour as a function of the 
centre-of-mass energy is in excellent agreement with the prediction. This gives confidence 
that this model can produce reliable predictions for the 6-hadron cross-section at the 
higher energies expected at the LHC. 

In Fig. [12] the cross-sections times dimuon branching fractions for the three T meson 



states are compared to the CSM NLO and NNLO* theoretical predictions 48 as a function 
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Figure 1 1 : Predictions based on the FONLL formalism 143 47 for the production cross-section 



for J/tp from b in the fiducial range < px < 14 GeV/c and 2.0 < y < 4.5 (yellow band). The 
uncertainty includes contributions from the renormalisation scale, quark masses and the choice 
of PDF set. The black dotted line shows the central value of the prediction. The points show 
the LHCb measurements at y/s = 2.76 7 [i~2], and 8 TeV. 



of pr- The NNLO* CSM provides a reasonable description of the experimental data, 
particularly for the T(3S) meson, which is expected to be less affected by feed-down. As 
for the prompt J/0 meson production, the CSM at NLO underestimates the cross-section 
by an order of magnitude. 



10 Conclusions 

The differential production cross-sections for J/0 and T mesons are measured as a function 
of pt an d y in the forward region, 2.0 < y < 4.5. The analysis is based on a data sample, 
corresponding to an integrated luminosity of 18 pb~ and 51 pb~ for the J/ip and T 
mesons, respectively, collected in the early part of 2012 at a centre-of-mass energy of 
a/s = 8 TeV. The production cross-sections of prompt J/0 mesons and J/0 from b are 
individually measured. An estimate of the bb total cross-section is also obtained. 

The results are compared with several recent theoretical predictions in the LHCb 
acceptance. The NNLO* CSM and the NLO NRQCD model (for the J/0 meson) provide 
a reasonable description of the experimental data on the production of prompt J/0 and 
T mesons as a function of under the assumption of zero polarisation. This confirms 
the progress in the theoretical calculations of quarkonium hadroproduction, even if the 
uncertainties on the predictions are still large. Theoretical predictions based on FONLL 
calculations are found to describe well the measured cross-section for J/0 from b and its 
dependence on centre-of-mass energy. 
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Figure 12: Comparison of the differential production cross-sections times dimuon branching 
fractions for (a) Y(1S), (b) T(2S) and (c) T(3S) mesons as a function of pt with direct production 
in an NNLO* CSM ^ (solid yellow) and an NLO CSM [7| (blue vertical shading) model. The 
points show the measurements reported in this analysis. 
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Table 4: Double-differential cross-section d PT dy m n b/(GeV/c) for prompt J/tp meson production 
in bins of of pp and y, with the assumption of no polarisation. The first error is statistical, the 
second is the component of the systematic uncertainty that is uncorrelated between bins and the 
third is the correlated component. 



MGeV/c) 



2.0 < y < 2.5 



2.5 < y < 3.0 



3.0 < y < 3.5 



0- 


1 


727.45 ± 13.14 ± 25.66 ± 


52.15 


772.40 ± 6.84 ± 11.94 ± 


55.38 


727.58±5.57±10.31±52.16 


1- 


2 


1463.31 ± 16.73 ± 19.26 ± 104.91 


1473.44 ± 8.74 ± 13.74 ± 105.64 


1357.44 ± 7.01 ± 11.10 ± 97.32 


2- 


3 


1237.49 ± 12.28 ± 22.39 ± 


88.72 


1196.41 ± 6.69 ± 20.45 ± 


85.78 


1083.40 ± 5.33 ± 


9.91 ±77.67 


3- 


4 


761.17 ± 


7.29 ± 20.60 ± 


54.57 


738.03 ± 4.20 ± 


8.94± 


52.91 


651.06 ± 3.45 ± 


6.62 ±46.68 


4- 


5 


432.67 ± 


4.33 ± 20.79 ± 


31.02 


408.62 ± 2.56 ± 


6.57± 


29.30 


359.82 ± 2.20 ± 


4.36 ±25.80 


5- 


6 


231.63 ± 


2.60 ± 


2.89 ± 


16.61 


217.04 ± 1.61 ± 


5.01 ± 


15.56 


183.07 ± 1.38 ± 


2.95 ±13.12 


6- 


7 


126.50 ± 


1.65 ± 


2.29 ± 


9.07 


116.27±1.07± 


4.31 ± 


8.34 


97.25±0.94± 


1.76 ± 6.97 


7- 


8 


68.05 ± 


1.07 ± 


1.98 ± 


4.88 


63.25 ± 0.74 ± 


1.09 ± 


4.53 


51.21 ± 0.64 ± 


1.23 ± 3.67 


8- 


9 


39.23 ± 


0.74 ± 


1.88 ± 


2.81 


34.85 ± 0.52 ± 


0.84± 


2.50 


27.34 ± 0.45 ± 


0.98± 1.96 


9- 


10 


22.04 ± 


0.52 ± 


0.66 ± 


1.58 


19.54 ± 0.37 ± 


0.49 ± 


1.40 


15.08 ± 0.32 ± 


0.70 ± 1.08 


10- 


11 


13.60 ± 


0.39 ± 


0.45 ± 


0.98 


11.59 ± 0.28 ± 


0.37± 


0.83 


8.88 ± 0.25 ± 


0.15± 0.64 


11- 


12 


8.06 ± 


0.28 ± 


0.23 ± 


0.58 


7.29 ± 0.22 ± 


0.33 ± 


0.52 


5.03±0.18± 


0.11 ± 0.36 


12- 


13 


5.26 ± 


0.22 ± 


0.20 ± 


0.38 


4.31 ± 0.16 ± 


0.06 ± 


0.31 


3.30±0.14± 


0.10± 0.24 


13- 


14 


3.30 ± 


0.17 ± 


0.21 ± 


0.24 


2.94±0.14± 


0.07± 


0.21 


2.09 ± 0.11 ± 


0.09± 0.15 






3.5 < y < 4.0 




4.0 < y 


< 4.5 








0- 


1 


670.50 ± 


5.34 ± 


6.78 ± 


48.07 


579.84 ± 6.08 ± 


7.65 ± 


41.57 






1- 


2 


1213.97 ± 


6.52 ± 20.90 ± 


87.04 


1003.19±7.18±10.37± 


71.92 






2- 


3 


926.23 ± 


4.88 ± 


7.50 ± 


66.41 


753.31 ± 5.75 ± 


6.80 ± 


54.01 






3- 


4 


542.09 ± 


3.19 ± 


5.56 ± 


38.87 


418.29 ± 3.89 ± 


4.61 ± 


29.99 






4- 


5 


284.82 ± 


2.00 ± 


3.44 ± 


20.42 


204.28 ± 2.31 ± 


2.68 ± 


14.65 






5- 


6 


146.41 ± 


1.29 ± 


2.20 ± 


10.50 


96.52 ± 1.39 ± 


2.61 ± 


6.92 






6- 


7 


73.68 ± 


0.84 ± 


1.33 ± 


5.28 


48.40 ± 0.90 ± 


0.92 ± 


3.47 






7- 


8 


37.39 ± 


0.56 ± 


0.82 ± 


2.68 


23.95 ± 0.59 ± 


0.60 ± 


1.72 






8- 


9 


20.05 ± 


0.40 ± 


0.52 ± 


1.44 


11.83 ± 0.39 ± 


0.34± 


0.85 






9- 


10 


11.04 ± 


0.29 ± 


0.38 ± 


0.79 


6.64±0.27± 


0.24± 


0.48 






10- 


11 


6.24 ± 


0.20 ± 


0.28 ± 


0.45 


3.40 ± 0.18 ± 


0.15 ± 


0.24 






11- 


12 


3.85 ± 


0.18 ± 


0.16± 


0.28 


2.05±0.14± 


0.10± 


0.15 






12- 


13 


2.23 ± 


0.13± 


0.05 ± 


0.16 


1.04 ± 0.09 ± 


0.03 ± 


0.07 






13- 


14 


1.49 ± 


0.10 ± 


0.04 ± 


0.11 


0.45 ± 0.06 ± 


0.02 ± 


0.03 
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Table 5: Double-differential cross-section d p T ^ in nb/(GeV/c) for the production of J/ip from b 
in bins of pr and y. The first error is statistical, the second is the component of the systematic 
uncertainty that is uncorrelated between bins and the third is the correlated component. 



p T (GeV/c) 


2.0 <y< 2.5 




2.5 < y < 3.0 




o.u < y < o.o 


0- 


1 


71.82 ± 2.53 ± 9.44 ± 


5.15 


71.70 ± 1.33 ± 6.58 ± 


5.14 


Ol.Oo zt l.lUzt0.0ozt4.4z 


1- 


2 


164.48 ±3.49 ±3.71 ±11.79 


157.59 ±1.84 ±3.88 ±11.30 


loo. zz zt i.4o zt z.oo zt y.oy 


2- 


3 


162.88 ±2.95 ±3.70 ±11.68 


152.06 ±1.62 ±3.38 ±10.90 


101 RQ _L 1 07-UOHS-US 70 
IZl.uo zt l.Z/ zt Z.Uo zt o. / Z 


3- 


4 


117.14±2.02±3.49± 


8.40 


106.89 ± 1.16 ± 1.79 ± 


7.66 


Q7 loin Q^X 1 KC\ 4- fi OK 
oi .10 zt u.yo zt 1.0U zt O.zO 


4- 


5 


75.00 ± 1.34 ± 3.71 ± 


5.38 


68.17±0.81±1.29± 


4.89 


coin C7_in 07 _|_ q oa 
Oo.uoztU.Of ztU.o/ zt O.04 


5- 


6 


46.32 ± 0.91 ± 0.97 ± 


3.32 


41.94±0.57±1.07± 


3.01 


qi a a i n ao _i_ a a a _l n OR 
ol.40 zt U.4o zt U.04 zt Z.ZO 


6- 


7 


28.96 ± 0.64 ± 0.56 ± 


2.08 


25.49 ± 0.42 ± 0.96 ± 


1.83 


1 o on _L A qr _l n Q7 -U 1 QI 
lo.oU zt U.oO zt U.O / zt l.ol 


7- 


8 


18.59 ± 0.46 ± 0.57 ± 


1.33 


15.36 ± 0.31 ± 0.31 ± 


1.10 


11 or,±n OA -U H 07 _l_ O. S21 
ll.ZO zt U.ZO zt U.Z / zt U.ol 


8- 


9 


11.40 ± 0.34 ± 0.55 ± 


0.82 


10.20±0.25±0.27± 


0.73 


£ QQ _|_ A OA I A o^-l-H AQ 

O.oo zt U.ZU zt U.ZO zt u.4y 


9- 


10 


7.77 ± 0.27 ± 0.24 ± 


0.56 


6.26±0.19±0.17± 


0.45 


4.04 zt U.lOzt U.Z1 zt U.OO 


10- 


11 


5.23±0.22±0.18± 


0.38 


4.16±0.15±0.13± 


0.30 


o sq xn io±n A^; _u a on 

Z.oo ztU.lozt U.UO zt U.ZU 


11- 


12 


3.46±0.17±0.10± 


0.25 


2.82±0.13±0.14± 


0.20 


1 7 c; _u a 1 a _l A A/1 -U n 1 Q 
l./OztU.lUztU. U4 zt U . 1 


12- 


13 


2.61±0.14±0.11± 


0.19 


2.17±0.11±0.03± 


0.16 


1 or,±n no -4- n n^ xn in 
l.oO zt u.uy zt U.UO zt U.1U 


13- 


14 


1.76±0.11±0.12± 


0.13 


1.39 ± 0.09 ± 0.04 ± 


0.10 


A OC J_ A A 7 _|_ A A/1 _|_ A flfi 

U.oOztU.Uf ztU.U4ztU.UO 






3.5 < y < 4.0 




4.0 < y < 4.5 






0- 


1 


46.51 ±1.04± 1.11 ± 


3.33 


29.15 ± 1.03 ± 1.41 ± 


2.09 




1- 


2 


100.11 ± 1.29 ± 3.76 ± 


7.18 


67.76 ± 1.35 ± 2.77 ± 


4.86 




2- 


3 


90.70 ± 1.20 ± 1.40 ± 


6.50 


58.35 ± 1.22 ± 0.72 ± 


4.18 




3- 


4 


62.80 ± 0.89 ± 1.38 ± 


4.50 


37.05 ± 0.90 ± 0.49 ± 


2.66 




4- 


5 


37.98 ± 0.59 ± 0.62 ± 


2.72 


21.13 ± 0.61 ± 0.44 ± 


1.52 




5- 


6 


20.61 ± 0.41 ±0.37± 


1.48 


12.04 ± 0.42 ± 0.34 ± 


0.86 




6- 


7 


12.09 ± 0.30 ± 0.22 ± 


0.87 


6.38 ± 0.29 ± 0.12 ± 


0.46 




7- 


8 


7.27±0.22±0.16± 


0.52 


3.69 ± 0.20 ± 0.09 ± 


0.26 




8- 


9 


4.19±0.17±0.12± 


0.30 


2.27±0.15±0.07± 


0.16 




9- 


10 


2.56±0.13±0.09± 


0.18 


1.09 ± 0.12 ± 0.09 ± 


0.08 




10- 


11 


1.42 ± 0.09 ± 0.07 ± 


0.10 


0.72±0.08±0.04± 


0.05 




11- 


12 


1.09 ± 0.08 ± 0.05 ± 


0.08 


0.37±0.05±0.04± 


0.03 




12- 


13 


0.73 ± 0.07 ± 0.02 ± 


0.05 


0.28 ± 0.04 ± 0.02 ± 


0.02 




13- 


14 


0.42 ± 0.05 ± 0.02 ± 


0.03 


0.16±0.03±0.01± 


0.01 
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Table 6: Fraction of J/ip from b (in %) in bins of of pr and y. The first uncertainty is statistical 
and the second systematic (uncorrected between bins). 



p T (GeV/c) 


on / „ ^ o k 
Z.O < y < z.o 


z.o < y < o.U 


o.U < y < o.o 


0- 


1 


8.9±0.3±1.1 


8.4±0.1±0.8 


7.7±0.1±0.8 


1- 


2 


10.1±0.2±0.2 


9.6 ±0.1 ±0.2 


9.0 ±0.1 ±0.2 


2- 


3 


11.6±0.2±0.2 


11.2±0.1±0.2 


10.0 ±0.1 ±0.1 


3- 


4 


13.3±0.2±0.2 


12.6 ±0.1 ±0.1 


11.7±0.1±0.2 


4- 


5 


14.7±0.2±0.2 


14.2 ±0.1 ±0.1 


12.9 ±0.1 ±0.1 


5- 


6 


16.6±0.3±0.3 


16.1 ±0.2 ±0.2 


14.5±0.2±0.2 


6- 


7 


1 o c i n q i n i 

lO.U —1— VJ . tj _1_ Vy . J- 


178J.nQJ.01 


1 C C J_ Q Q 1 Q 


7- 


8 


21.3±0.4±0.2 


19.3±0.3±0.2 


17.7±0.3±0.1 


8- 


9 


22.4±0.5±0.1 


22.4±0.5±0.2 


19.7±0.5±0.1 


9- 


10 


25.8±0.7±0.2 


23.9±0.6±0.2 


22.5±0.7±0.2 


10- 


11 


27.3±0.9±0.2 


25.9±0.8±0.1 


23.4±0.9±0.3 


11- 


12 


29.6±1.1±0.2 


27.2 ±1.0 ±0.5 


24.9 ±1.2 ±0.1 


12- 


13 


32.3 ±1.3 ±0.6 


32.8 ±1.3 ±0.2 


27.8 ±1.5 ±0.5 


13- 


14 


34.1±1.7± 1.0 


30.9 ±1.6 ±0.4 


26.9±2.0±0.4 






o.o < y < 4.u 


4.U < y < 4.0 




0- 


1 


6.4±0.1±0.1 


4.6±0.2±0.2 




1- 


2 


7.5±0.1±0.3 


6.2±0.1±0.2 




2- 


3 


8.8±0.1±0.1 


7.1 ±0.1 ±0.1 




3- 


4 


10.3±0.1±0.2 


8.0 ±0.2 ±0.1 




4— 


5 


11.6±0.2±0.1 


9.3 ±0.2 ±0.1 




5- 


6 


12.1±0.2±0.1 


10.9±0.4±0.1 




6- 


7 


13.8±0.3±0.1 


11.4±0.5±0.0 




7- 


8 


16.0±0.4±0.1 


13.2±0.7±0.1 




8- 


9 


16.8±0.6±0.2 


15.8±0.9±0.2 




9- 


10 


18.3±0.8±0.1 


13.2 ±1.3 ±0.9 




10- 


11 


18.1 ±1.1 ±0.3 


16.7±1.7±0.6 




11- 


12 


21.5±1.5±0.2 


14.9 ±2.0 ±1.3 




12- 


13 


23.3±2.0±0.4 


20.3 ±3.0 ±1.1 




13- 


14 


21.2±2.1±0.6 


25.8±4.4± 1.9 
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Table 7: Double-differential production cross-sections d ^ T j x B IS in pb/(GeV/c) for the T(1S) 
meson in bins of transverse momentum and rapidity, assuming no polarisation. The first error 
is statistical, the second is the component of the systematic uncertainty that is uncorrelated 
between bins and the third is the correlated component. 



MGeV/c) 


2.0 <y< 2.5 


2.5 < y < 3.0 


3.0 < y < 3.5 


0- 


-1 


73.9 ± 4.2 ± 1.8 ± 5.3 


64.7 ± 2.6 ± 1.3 ± 4.6 


62.1 ± 2.6 ± 0.7 ± 4.4 


1- 


-2 


197.7 ± 6.8 ± 2.7 ± 14.1 


185.6 ± 4.5 ± 1.9 ± 13.2 


157.2 ± 4.2 ± 1.7 ± 11.2 


2- 


-3 


248.1 ± 7.5 ± 2.7 ± 17.6 


228.5 ± 5.0 ± 2.1 ± 16.2 


210.6 ± 4.9 ± 1.6 ± 15.0 


3 


-4 


237.9 ± 7.2 ± 4.1 ± 16.9 


243.5 ± 5.1 ± 1.5 ± 17.3 


210.3 ± 4.9 ± 2.2 ± 15.0 


4- 


-5 


220.6 ± 6.9 ± 3.1 ± 15.7 


199.1 ± 4.7 ± 1.2 ± 14.2 


185.0 ± 4.6 ± 0.8 ± 13.2 


5- 


-6 


181.2 ± 6.1 ± 4.2 ± 12.9 


177.4 ± 4.4 ± 1.8 ± 12.6 


161.4 ± 4.3 ± 0.8 ± 11.5 


6- 


-7 


140.1 ± 5.3 ± 4.5 ± 10.0 


128.2 ± 3.7 ± 0.9 ± 9.1 


116.0 ± 3.7 ± 0.6 ± 8.3 


7- 


-8 


106.2 ± 4.6 ± 1.6 ± 7.6 


106.3 ± 3.4 ± 1.5 ± 7.6 


93.9 ± 3.3 ± 1.2 ± 6.7 


8- 


-9 


75.3 ± 3.8 ± 2.3 ± 5.4 


74.0 ± 2.8 ± 0.7 ± 5.3 


62.3 ± 2.6 ± 0.5 ± 4.4 


9- 


-10 


60.9 ± 3.3 ± 1.5 ± 4.3 


56.7 ± 2.4 ± 0.7 ± 4.0 


51.1 ± 2.4 ± 0.4 ± 3.6 


10- 


-11 


51.4 ± 3.0 ± 1.7 ± 3.7 


44.4 ± 2.1 ± 0.8 ± 3.2 


32.9 ± 1.9 ± 0.3 ± 2.3 


11- 


-12 


34.7 ± 2.5 ± 1.4 ± 2.5 


31.1 ± 1.7 ± 0.6 ± 2.2 


25.6 ± 1.6 ± 0.3 ± 1.8 


12- 


-13 


29.7 ± 2.2 ± 0.3 ± 2.1 


21.5 ± 1.4 ± 0.9 ± 1.5 


17.0 ± 1.3 ± 0.2 ± 1.2 


13 


-14 


22.3 ± 1.9 ± 0.5 ± 1.6 


16.8 ± 1.3 ± 0.1 ± 1.2 


12.7 ± 1.1 ± 0.2 ± 0.9 


14- 


-15 


15.6 ± 1.6 ± 0.8 ± 1.1 


13.8 ± 1.1 ± 0.5 ± 1.0 


10.4 ± 1.0 ± 0.2 ± 0.7 






3.5 < y < 4.0 


4.0 < y < 4.5 




0- 


-1 


52.0 ± 3.2 ± 0.8 ± 3.7 


39.7 ± 5.8 ± 0.4 ± 2.8 




1 

1- 


— z 


128.0 ± 4.9 ± 0.9 ± 9.1 


82.6 ± 8.0 ± 1.9 ± 5.9 




2- 


-3 


167.1 ± 5.5 ± 2.2 ± 11.9 


125.4 ± 9.6 ± 1.3 ± 8.9 




3 


-4 


166.4 ± 5.5 ± 1.2 ± 11.8 


130.2 ± 9.3 ± 1.2 ± 9.3 




4- 


-5 


148.9 ± 5.2 ± 1.3 ± 10.6 


95.6 ± 7.6 ± 1.1 ± 6.8 




5- 


-6 


101.2 ± 4.2 ± 0.5 ± 7.2 


68.4 ± 6.2 ± 0.6 ± 4.9 




6- 


-7 


79.1 ± 3.7 ± 0.5 ± 5.6 


58.5 ± 5.6 ± 0.5 ± 4.2 




7- 


-8 


63.9 ± 3.3 ± 0.6 ± 4.5 


41.2 ± 4.7 ± 0.5 ± 2.9 




8- 


-9 


47.0 ± 2.8 ± 0.3 ± 3.3 


30.5 ± 3.9 ± 0.6 ± 2.2 




9- 


-10 


34.8 ± 2.4 ± 0.3 ± 2.5 


16.5 ± 2.8 ± 0.2 ± 1.2 




10- 


-11 


23.5 ± 2.0 ± 0.2 ± 1.7 


13.4 ± 2.4 ± 0.2 ± 1.0 




11- 


-12 


17.8 ± 1.7 ± 0.2 ± 1.3 


12.1 ± 2.3 ± 0.2 ± 0.9 




12- 


-13 


13.4 ± 1.4 ± 0.2 ± 1.0 


6.6 ± 1.7 ± 0.1 ± 0.5 




13 


-14 


10.0 ± 1.2 ± 0.2 ± 0.7 


3.4 ± 1.1 ± 0.1 ± 0.2 




14- 


-15 


7.2 ± 1.0 ± 0.1 ± 0.5 


0.7 ± 0.5 ± 0.0 ± 0.1 
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Table 8: Double-differential production cross-sections d p T d y x $ 2S m pb/(GeV/c) for the T(2S) 
meson in bins of transverse momentum and rapidity, assuming no polarisation. The first error 
is statistical, the second is the component of the systematic uncertainty that is uncorrelated 
between bins and the third is the correlated component. Regions where the number of events 
was not large enough to perform a measurement are indicated with a dash. 



p T ( GeV/c) 2.0 < y < 2.5 2.5 < y < 3.0 3.0 < y < 3.5 



0- 


-1 


14.5 


± 


1.8 


± 


0.3 


± 


1.0 


15.2 


± 


1.3 


± 


0.7 


± 


1.1 


11.7 


± 


1.1 


± 


0.1 


± 


0.8 


1- 


-2 


33.2 


± 


2.7 


± 


1.3 


± 


2.4 


42.6 


± 


2.1 


± 


0.3 


± 


3.1 


39.0 


± 


2.1 


± 


0.2 


± 


2.8 


2- 


-3 


49.3 


± 


3.3 


± 


1.0 


± 


3.5 


54.3 


± 


2.4 


± 


1.0 


± 


3.9 


45.0 


± 


2.2 


± 


0.2 


± 


3.2 


3- 


-4 


56.8 


± 


3.5 


± 


0.6 


± 


4.1 


61.0 


± 


2.6 


± 


0.4 


± 


4.4 


48.2 


± 


2.3 


± 


0.2 


± 


3.5 


4- 


-5 


46.4 


± 


3.1 


± 


0.5 


± 


3.3 


45.1 


± 


2.2 


± 


0.5 


± 


3.2 


39.5 


± 


2.1 


± 


0.3 


± 


2.8 


5- 


-6 


51.1 


± 


3.2 


± 


0.8 


± 


3.7 


48.9 


± 


2.3 


± 


0.3 


± 


3.5 


38.1 


± 


2.1 


± 


0.6 


± 


2.7 


6- 


-7 


43.4 


± 


2.9 


± 


0.9 


± 


3.1 


34.9 


± 


1.9 


± 


0.6 


± 


2.5 


29.3 


± 


1.8 


± 


0.2 


± 


2.1 


7- 


-8 


28.6 


± 


2.4 


± 


0.6 


± 


2.1 


29.3 


± 


1.8 


± 


0.7 


± 


2.1 


25.5 


± 


1.7 


± 


0.3 


± 


1.8 


8- 


-9 


21.2 


± 


2.0 


± 


0.1 


± 


1.5 


23.7 


± 


1.6 


± 


0.2 


± 


1.7 


17.3 


± 


1.4 


± 


0.1 


± 


1.2 


9- 


10 


20.9 


± 


2.0 


± 


0.2 


± 


1.5 


16.1 


± 


1.3 


± 


0.2 


± 


1.2 


13.2 


± 


1.2 


± 


0.1 


± 


0.9 


10- 


-11 


16.0 


± 


1.7 


± 


0.2 


± 


1.1 


13.4 


± 


1.2 


± 


0.2 


± 


1.0 


13.4 


± 


1.2 


± 


0.2 


± 


1.0 


11- 


-12 


14.9 


± 


1.6 


± 


0.2 


± 


1.1 


10.2 


± 


1.0 


± 


0.2 


± 


0.7 


10.1 


± 


1.0 


± 


0.1 


± 


0.7 



12- 


-13 


7.8 ± 1.2 ± 0.2 ± 0.6 


10.4 ± 1.0 ± 0.2 ± 0.7 


8.4 ± 1.0 ± 0.2 ± 0.6 


13- 


-14 


5.6 ± 1.0 ± 0.2 ± 0.4 


7.0 ± 0.8 ± 0.1 ± 0.5 


5.4 ± 0.8 ± 0.1 ± 0.4 


14- 


-15 


4.0 ± 0.8 ± 0.3 ± 0.3 


5.8 ± 0.8 ± 0.1 ± 0.4 


3.5 ± 0.6 ± 0.1 ± 0.3 






3.5 < y < 4.0 


4.0 < y < 4.5 




0- 


-1 


11.2 ± 1.5 ± 0.1 ± 0.8 






1- 


-2 


26.3 ± 2.2 ± 0.2 ± 1.9 






2- 


-3 


35.0 ± 2.5 ± 0.3 ± 2.5 






3- 


-4 


44.7 ± 2.9 ± 0.3 ± 3.2 






4- 


-5 


40.6 ± 2.7 ± 0.2 ± 2.9 


14.5 ± 2.9 ± 0.1 ± 1.0 




5- 


-6 


28.4 ± 2.3 ± 0.1 ± 2.0 


12.8 ± 2.8 ± 0.1 ± 0.9 




6- 


-7 


19.8 ± 1.8 ± 0.1 ± 1.4 


19.5 ± 3.2 ± 0.3 ± 1.4 




7- 


-8 


20.1 ± 1.9 ± 0.1 ± 1.4 


17.0 ± 3.0 ± 0.2 ± 1.2 




8- 


-9 


14.8 ± 1.6 ± 0.3 ± 1.1 


8.7 ± 2.1 ± 0.1 ± 0.6 




9- 


10 


8.7 ± 1.2 ± 0.1 ± 0.6 






10- 


-11 


6.5 ± 1.1 ± 0.2 ± 0.5 






11- 


-12 


7.0 ± 1.1 ± 0.1 ± 0.5 






12- 


-13 


2.7 ± 0.7 ± 0.1 ± 0.2 






13- 


-14 


4.3 ± 0.8 ± 0.1 ± 0.3 






14- 


-15 
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Table 9: Double-differential production cross-sections d p T d y x ^ S m pb/(GeV/c) for the T(35) 
meson in bins of transverse momentum and rapidity, assuming no polarisation. The first error 
is statistical, the second is the component of the systematic uncertainty that is uncorrelated 
between bins and the third is the correlated component. Regions where the number of events 
was not large enough to perform a measurement are indicated with a dash. 



MGeV/c) 


z.u <~. y <. z.o 


Z.O <^ y <, O.U 


O.U o.O 


0- 


-1 


4.8 ± 1.0 ± 0.1 ± 0.4 


5.5 ± 0.8 ± 0.1 ± 0.4 


4.7 ± 0.7 ± 0.1 ± 0.3 


1- 


-2 


17.9 ± 2.0 ± 0.3 ± 1.3 


18.2 ± 1.4 ± 0.2 ± 1.3 


17.3 ± 1.4 ± 0.1 ± 1.3 


2- 


-3 


25.5 ± 2.3 ± 0.4 ± 1.9 


17.3 ± 1.3 ± 0.1 ± 1.3 


21.3 ± 1.5 ± 0.1 ± 1.6 


3- 


-4 


22.9 ± 2.2 ± 0.3 ± 1.7 


26.2 ± 1.7 ± 0.2 ± 1.9 


22.8 ± 1.6 ± 0.1 ± 1.7 


4- 


-5 


23.6 ± 2.2 ± 0.3 ± 1.7 


18.0 ± 1.4 ± 0.1 ± 1.3 


19.9 ± 1.5 ± 0.1 ± 1.5 


5- 


-6 


19.6 ± 2.0 ± 0.3 ± 1.4 


21.8 ± 1.5 ± 0.2 ± 1.6 


18.4 ± 1.4 ± 0.1 ± 1.4 


6- 


-7 


21.4 ± 2.1 ± 0.4 ± 1.6 


18.9 ± 1.4 ± 0.2 ± 1.4 


19.8 ± 1.5 ± 0.2 ± 1.5 


7- 


-8 


13.7 ± 1.6 ± 0.3 ± 1.0 


16.1 ± 1.3 ± 0.2 ± 1.2 


12.0 ± 1.2 ± 0.2 ± 0.9 


8- 


-9 


14.2 ± 1.6 ± 0.6 ± 1.1 


12.5 ± 1.1 ± 0.2 ± 0.9 


9.3 ± 1.0 ± 0.1 ± 0.7 


9- 


-10 


8.7 ± 1.3 ± 0.3 ± 0.6 


10.4 ± 1.0 ± 0.2 ± 0.8 


10.8 ± 1.1 ± 0.1 ± 0.8 


10- 


-11 


9.6 ± 1.3 ± 0.5 ± 0.7 


7.4 ± 0.9 ± 0.3 ± 0.5 


6.8 ± 0.8 ± 0.1 ± 0.5 


11- 


-12 


8.1 ± 1.2 ± 0.1 ± 0.6 


5.9 ± 0.8 ± 0.2 ± 0.4 


4.0 ± 0.7 ± 0.1 ± 0.3 


12- 


-13 


6.3 ± 1.0 ± 0.1 ± 0.5 


5.3 ± 0.7 ± 0.1 ± 0.4 


3.0 ± 0.6 ± 0.0 ± 0.2 


13- 


-14 


10.0 ± 1.3 ± 0.2 ± 0.7 


4.6 ± 0.7 ± 0.2 ± 0.3 


3.7 ± 0.6 ± 0.1 ± 0.3 


14- 


-15 


4.0 ± 0.8 ± 0.1 ± 0.3 


2.2 ± 0.5 ± 0.1 ± 0.2 


3.5 ± 0.6 ± 0.1 ± 0.3 






3.5 < y < 4.0 


4.0 < y < 4.5 




0- 


-1 


3.3 ± 0.8 ± 0.0 ± 0.2 


- 




1- 


-2 


19.1 ± 1.9 ± 0.2 ± 1.4 


- 




2- 


-3 


20.0 ± 1.9 ± 0.2 ± 1.5 






3- 


-4 


18.6 ± 1.8 ± 0.1 ± 1.4 






4- 


-5 


18.5 ± 1.8 ± 0.1 ± 1.4 


6.3 ± 2.0 ± 0.1 ± 0.5 




5- 


-6 


11.0 ± 1.4 ± 0.1 ± 0.8 


5.3 ± 1.7 ± 0.1 ± 0.4 




6- 


-7 


9.3 ± 1.3 ± 0.1 ± 0.7 


1.0 ± 0.7 ± 0.0 ± 0.1 




7- 


-8 


8.1 ± 1.2 ± 0.1 ± 0.6 


7.3 ± 1.9 ± 0.1 ± 0.5 




8- 


-9 


8.0 ± 1.2 ± 0.1 ± 0.6 


3.0 ± 1.2 ± 0.0 ± 0.2 




9- 


-10 


7.9 ± 1.1 ± 0.1 ± 0.6 






10- 


-11 


3.6 ± 0.8 ± 0.0 ± 0.3 






11- 


-12 


5.1 ± 0.9 ± 0.1 ± 0.4 






12- 


-13 


4.5 ± 0.9 ± 0.1 ± 0.3 






13- 


-14 


0.7 ± 0.3 ± 0.0 ± 0.1 






14- 


-15 
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Table 10: Ratios of cross-sections T(2S) — > fi + fx~ and T(3S) — > ^ + ^~ with respect to 
T(1S) — > as a function of in the range 2.0 < y < 4.0, assuming no polarisation. 

The first error is statistical, the second is the component of the systematic uncertainty that is 
uncorrelated between bins and the third is the correlated component. 



PT ( GeV/c) R 2S ' IS R 3S / 1S 



0- 


-1 


0. 


.210 


± 


0. 


.013 


± 





.004 


± 


0. 


.006 





.075 


± 


0. 


.007 


± 


0. 


.001 


± 


0, 


,003 


1- 


-2 


0. 


.219 


± 


0. 


.008 


± 





.002 


± 





.006 


0. 


.109 


± 


0. 


.005 


± 


0. 


.001 


± 


0, 


,004 


2- 


-3 





.219 


± 


0. 


.007 


± 





.002 


± 





.006 


0. 


.096 


± 


0. 


.004 


± 


0. 


,001 


± 


0, 


003 


3- 


-4 


0, 


,244 


± 


0. 


.008 


± 


0. 


.002 


± 





.007 





.107 


± 


0. 


.004 


± 





,001 


± 


0, 


,004 


4- 


-5 


0. 


,226 


± 


0. 


,008 


± 


0. 


.001 


± 


0. 


.006 


0. 


.104 


± 


0. 


.005 


± 


0. 


,001 


± 


0, 


,004 


5- 


-6 


0. 


,267 


± 


0. 


,009 


± 


0. 


.003 


± 


0. 


.007 


0. 


.116 


± 


0. 


.005 


± 





,001 


± 


0, 


,004 


6- 


-7 


0. 


,270 


± 


0. 


,011 


± 


0. 


.004 


± 


0. 


.008 


0. 


.151 


± 


0. 


.007 


± 


0. 


,001 


± 


0, 


,005 


7- 


-8 





,280 


± 


0. 


,012 


± 


0. 


.003 


± 


0. 


.008 


0. 


.137 


± 


0. 


.007 


± 


0. 


,001 


± 


0, 


,005 


8- 


-9 


0. 


,298 


± 


0. 


,015 


± 





.003 


± 


0. 


.008 


0. 


.170 


± 


0. 


010 


± 


0. 


,002 


± 


0, 


,006 


9- 


10 


0. 


,280 


± 


0. 


,016 


± 


0. 


.003 


± 


0. 


.008 





.190 


± 


0. 


.012 


± 


0. 


,002 


± 


0, 


,007 


10- 


-11 





,330 


± 


0. 


,021 


± 


0. 


.005 


± 


0. 


.009 


0. 


.181 


± 


0. 


013 


± 





,004 


± 


0, 


,006 


11- 


-12 


0. 


,379 


± 


0. 


,027 


± 





.006 


± 





.011 





.203 


± 


0. 


.017 


± 


0. 


,003 


± 


0, 


,007 


12- 


-13 





,401 


± 


0. 


,032 


± 





.007 


± 





.011 


0. 


.235 


± 


0. 


.021 


± 


0. 


003 


± 


0, 


,008 


13- 


-14 


0. 


.388 


± 


0. 


.036 


± 


0. 


.006 


± 


0. 


.011 


0. 


.292 


± 


0. 


.027 


± 


0. 


,004 


± 


0, 


,010 



Table 11: Ratios of cross-sections Y(2S) — > ^ + ^~ and T(3S) — > ^ + ^~ with respect to 
T(1S) — > fi + as a function of y in the range 2.0 < y < 4.0, assuming no polarisation. The first 
error is statistical, the second is the component of the systematic uncertainty that is uncorrelated 
between bins and the third is the correlated component. 



y 



p>2S/lS fi?,S/li 



2.0-2.5 0.2545 ± 0.0068 ± 0.0021 ± 0.0071 0.1260 ± 0.0041 ± 0.0008 ± 0.0044 

2.5-3.0 0.2653 ± 0.0050 ± 0.0014 ± 0.0075 0.1210 ± 0.0029 ± 0.0005 ± 0.0043 

3.0-3.5 0.2476 ± 0.0051 ± 0.0009 ± 0.0070 0.1251 ± 0.0032 ± 0.0004 ± 0.0044 

3.5-4.0 0.2558 ± 0.0075 ± 0.0011 ± 0.0072 0.1315 ± 0.0048 ± 0.0004 ± 0.0046 
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